Statistical mechanics of macromolecular networks without rephcas 



oo 



O 
O 



o 



Michael P Solf, Thomas A Vilgis 
Max- Planck- Institut fur Polymerforschung, Postfach 3148, 55021 Mainz, Germany 
(to appear in J. Phys. A: Math. Gen. 28, December (1995)) 

Abstract 

We report on a novel approach to the Deam-Edwards model for interacting 
polymeric networks without using replicas. Our approach utilizes the fact that 
a network modelled from a single non-interacting Gaussian chain of macro- 



> 

\Q ' scopic size can be solved exactly, even for randomly distributed crosslinking 



junctions. We derive an exact expression for the partition function of such a 
generalized Gaussian structure in the presence of random external fields and 



' for its scattering function 5*0. We show that 5*0 of a randomly crosslinked 

a" 

Gaussian network (RCGN) is a self-averaging quantity and depends only 



on crosslink concentration M /N, where M and N are the total numbers of 
crosslinks and monomers. From our derivation we find that the radius of gy- 



^ ■ ration of a RCGN is of the universal form = (0.26 ± 0.01)a^N/M, with 



a being the Kuhn length. To treat the excluded volume effect in a system- 
atic, perturbative manner, we expand the Deam-Edwards partition function 
in terms of density fiuctuations analogous to the theory of linear polymers. 
For a highly crosslinked interacting network we derive an expression for the 
free energy of the system in terms of 5*0 which has the same role in our model 
as the Debye function for linear polymers. Our ideas are easily generalized 
to crosslinked polymer blends which are treated within a modified version of 
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Leibler's mean field theory for block copolymers. 
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I. INTRODUCTION 



Randomly crosslinked macromolecules present a challenging field from a physical as well 
as from a mathematical point of view with many practical applications in polymer sciences. 
For this it is very unfortunate that the statistical mechanics of polymer networks is still 
poorly understood. Although there has been a great deal of theoretical interest in this topic 
over the last couple of years we feel that a satisfying answer of how to deal with 

polymer networks from a microscopic point of view is still missing. One of the reasons is 
that most of the recent work on the subject is based on replica field theory originally 
introduced into polymer science by Edwards et al in which, however, much of the 
underlying physics remains hidden in the complicated mathematics of the replica formalism 
i- 

It is the purpose of this paper to present an alternative approach to the statistical 
mechanics of randomly crosslinked macromolecules that goes beyond the earlier phantom 
type models (for details see, for example, reference but at the same time avoids the well 
known difficulties associated with the replica trick [^]. Our theory is based on the minimal 
network model by Deam and Edwards however, avoids replica field theory completely. 

The strategy in this paper is outlined below: As generalization of the Wiener measure 
in the theory of linear polymers []TT| we introduce the concept of a randomly crosslinked 
Gaussian network (RCGN) in section 2. As we will show in section 3, the non-interacting 
problem can be solved exactly, even for random crosslinking junctions. In section 3 we 
develop our general formalism and derive an exact formula for the partition function of a 
generalized Gaussian structure in the presence of random external fields. This equation 
(20) is the central mathematical result of the paper, for it is also a generating function 
from which further results are obtained. In section 4 we consider some applications. In 
particular we calculate the static structure function 5*0 for RCGNs without excluded volume 
interaction. The important finding here is that 5*0 is a self- aver aging quantity, i.e., it does not 
depend on the topological details of the model. The consequence is that for RCGNs 5*0 is a 
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quasi universal function that has a similar role for polymer networks as the Debye function 
for linear polymers. The interacting case is treated in section 5. To take into account 
excluded volume interaction, we transform our original network Hamiltonian to collective 
density variables. We work here in close analogy with the excluded volume problem for 
dense polymer melts [TT[]. For highly crosslinked systems it is sufficient to consider only 



lowest order density fluctuations, although higher order terms are readily calculated within 
our formalism. In this case it is easy to show that the free energy F can be expressed in 
terms of the scattering function 5*0 of the non-interacting system and the excluded volume 
parameter. It is crucial to realize that F depends on the frozen degrees of freedom (random 
crosslinking junctions) only via 5*0, and no further quenched averaging remains to be done. 
To treat crosslinked polymer blends we apply a modifled version of Leibler's theory for 
microphase separation in block copolymers |12[ and show how to obtain similar criteria for 
phase instability in mult i- component networks. Generalizations and outlook are given in 
section 6. However, more detailed calculations of the latter applications will be presented 
in a forthcoming publication. 

II. FORMULATION OF THE MODEL 

We consider flexible interacting macromolecules on the level of the Edwards Hamiltonian 
For a single polymer chain in d spatial dimensions the Edwards Hamiltonian consists of 
two parts H = + H^, where 

d ^ 

f3H^= (1) 

is the Wiener measure that models the connectedness of the chain as a Gaussian. Self- 
avoidance between monomers is described by a pseudopotential of the form 

N 

PH, = v <^(R.-R.) , (2) 

0<i<j 

where 5(R) is the Dirac delta function. In (1) and (2) we have adopted the following notation: 
f > is the second virial coefficient that characterizes the repulsive excluded volume 



interaction between monomers, is the degree of polymerization, a the Kuhn length, and 
Rj {i = 0, .., A^) are monomer coordinates; (3 = (kj^T)^^ as usual. For convenience we restrict 
the following discussion to networks that are modelled from one single but huge polymer 
chain of macroscopic dimension (figure 1). It was shown that this simplification gives the 
correct physics for highly crosslinked polymer networks above the percolation threshold [|T|. 
A generalization of our method to multi-polymer networks will be discussed later on in 
section 6. 

To describe M permanently crosslinked monomers, we specify each junction by a pair 
of randomly chosen "crosslink coordinates" ie,je (0 < ie,je ^ N, e = 1,...,M), such that 
monomer Rj^ is connected to monomer R^^ (figure 1). The whole set of junction points 
C = {ie,je}eLi represents the random connectivity of the network. Within the framework 
of the Deam-Edwards model the partition function of a Gaussian phantom network with 
excluded volume is given by 

„ N AI 

ZiC) = JY[dR, e-^(^w+/^i) -Q s (R^^ - R^.J , (3) 

Y i=0 e=l 

where total phase space is now restricted by the additional crosslinking constraints. Equation 
(3) describes a phantom network in a sense that the polymer chain is free to pass through 
itself irrespective of entanglements. Chain motion is only restricted by the presence of 
permanent crosslinks and the excluded volume interaction. 

In the replica formalism the next step is to perform a quenched average over the loga- 
rithm of Z{C) with a suitable distribution for the "frozen" crosslink coordinates {ie,je}eLi- 
Using the replica trick ||^ this leads to a non-trivial modification of the interaction term 
Hi in which all replicas become coupled Applying standard techniques for setting up 

field theories in polymer physics, we were able to map equation (3) into a dn dimensional 
0(m) field theory in the limit n,m —>■ Unfortunately in the replica formalism further 
progress highly relies upon crude approximations or variational assumptions 

Contrary to the replica method we do not carry out the quenched average at this stage 
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of the calculation. Instead we model the delta function in (3) by a Gaussian distribution 
with width e in the limit £ ^ 0. Therefore we are keeping all random crosslink coordinates 
explicitly in the partition function. Formally this means that if the Wiener measure (3H^ is 
replaced by the more general expression for a RCGN 

d ^ d ^ 

PHo = E(R^ - R-i)' + ^ E(R^e - R>)^ , (4) 

i=l e=l 

we can eliminate the delta constraint in (3). With equations (2) and (4) we are now in a 
position to introduce our RCGN Hamiltonian as follows 

H = Ho + H,. (5) 

To enforce the crosslinking constraints we finally have to perform the limit e — > 0. It is 
mathematically convenient to normalize the Gaussian measure defined by (4) with respect 
to the non-interacting system which will serve us as a reference state. In this case the 
partition function in (3) can be cast into the more convenient form 

ZiO = (e-^^^)^ = exp ( - PiF - Fo)) , (6) 

where the average (...)q stands for 

. N N 

Jim (/ n rfR. x{Jl[dR, e-^^«)" . (7) 

i=0 i=0 

For V = the free energy F of the interacting network reduces to that of a non-interacting 
RCGN denoted by Fq. The partition function (6) is completely equivalent to the one used 
by Edwards et al . It is also identical to the one used in more recent works by Goldbart 
and coworkers 0,^, except for the fact that in this paper we are dealing with a single-chain 
network exclusively. This is primarily a matter of convenience which has been discussed 
in greater detail in references The Hamiltonian (5) is of course easily generalized to 

multi-chain or multi-component networks. 

To make further analytic progress, it will become useful to decouple the interaction term 
Hi by rewriting (2) in terms of collective density variables $k = Y^f=o exp(2kRj) and applying 
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the Hubbard-Stratonovich transformation to (6). Since this is standard procedure @JTT| we 
quote here only the final result 



Z{C) oc/nk>o^^'fkexp(-^Ek>o'^>k$-kj (8) 
exp (z^E.=oi:k^o'^'-kexp(zkR,)))^ . 



X 

By k > we mean the positive half of k space, i.e., > 0, and V is the volume. The 
partition function (8) in combination with the measure defined by (7) is equivalent to the 
minimal network model of Deam and Edwards and will be our starting point for further 
investigations. 



III. GENERAL FORMALISM 

Before dealing with the more complicated excluded volume situation directly, it is easier 
to develop the mathematical formalism for non- interacting RCGNs at first. To evaluate 
equation (8) in a systematic fashion, it is standard to expand the first exponential in the 
entropic, second part of (8) in terms of density fiuctuations $_k or the excluded volume 
parameter. In any case the main mathematical task is to calculate averages which are of the 
general form ( exp(ib*R))Q, where we have introduced dx (A^+1) dimensional "supervectors" 
R = (Ro, Ri, Rtv)* and b = (bo, bi, bAr)*. For the moment we assume that b is 
completely arbitrary, but does not depend on R explicitly. By bjRj we will always mean a d 
dimensional inner vector product, and t denotes the transposed vector. Since (.■■)q involves 
a Gaussian integration with the measure defined by (7), it is essential to find the inverse 
matrix of the quadratic form, equation (4), first. For this it is most convenient to switch to 
matrix notation. 

Using matrix notation in which R is the above defined "supervector" , the quadratic form 
in (4) can be written as follows 

d ^'^ 

PHa = — R\zW + 5:/C(2e, Je))R , (9) 

e=l 
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where z — {s/af, and symmetric (A^ + 1) x (A^ + 1) matrices W and /C. Here 



^ 1-1 ••• 0^ 



-1 2 -1 ••• 

... -1 2-1 
••• 0-1 1 



V 



(10) 



denotes the connectivity matrix associated with i^w, and 



^0 0^ 



1 ; -1 
••• ••• 

-1 : 10 





ig-th position 



je-th position 



(11) 



represents a single crosshnking junction at random position {ie,je)- Equation (9) is easily 
verified by inspection. 

For the following it is essential to note that the quadratic form, equation (9), is only 
semzpositive definite, hence no matrix inverse does exist. The non- negativity of He follows 
immediately from equation (4), since it is a sum of squares. However, there is one zero 
eigenvalue associated with eigenvector R = (1, 1)*. This is easily seen from noting that 
the column sum of the matrix in (9) is always zero irrespective of the value of z. Before 
we proceed, we need to transform (9) to a positive definite quadratic form by removing 
the mode which belongs to eigenvalue zero. This mode corresponds to a displacement of 
the center of mass, and since we are dealing with a simply-connected structure, there can 
be only one such mode. The desired transformation is easily accomplished by switching to 
internal coordinates = Rj — R^-i, or in matrix notation R = T>r, where 
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^ Ro ^ 



r = 



and V — 



1 ••• 
1 1 ••• 

1 1 •■ 



(12) 



For the multi-chain network we refer to our discussion in section 6. With the above manip- 
ulations the measure corresponding to equation (9) takes on the simple form 

N , N 



^o(C) = / n ^R-i = V / n rfr, exp ( - — v'M r) , 

i=Q •' i=l 



(13) 



where r = (ri, rjv)*, and we performed an integration over Rq. In the limit £ — > equation 
(13) yields the partition function of a RCGN without excluded volume. The dimensional 
matrix A4 in (13) is given by 



M{z)^zI + VV' , 



(14) 



where I denotes an dimensional unit matrix. The outer matrix product VV^ is formed 
with the N X M "crosshnk matrix" V — (pi, Pm), where each column vector is defined 

by 



(0,...,0,1,1,....L1,0,...,0)* , e=l,...,M. 



(15) 



ie + l to je 

The I's in (15) run from the (ie -l- l)th to the jeth position, the rest of the components is 0. 
For convenience we will assume that jg is always larger than i^. Equation (14) is derived in 
appendix A. 

The crosslink matrix V is defined in such a way that it contains complete information 
about the crosslink positions in an unique way. By construction it is exactly this extra 
term in A4 that distinguishes the network problem from a linear polymer. For the following 
derivation it is crucial to have the crosslinking constraint in equation (14) in form of an 
outer matrix product of V. Note that A4 is positive definite, and thus its inverse exists. 

Using equations (13) and (14) the measure in (7) can be redefined in terms of M. 
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-I . . N J 



G?rj exp - - 

i=l 

In the remainder of this section we will show that this average exists for e —>■ and derive 
an exact expression. 

With (16) we are in a position to calculate expectation values of the form (exp(ib*R))Q. 
Since r depends only linearly on R, it is sufficient to consider expressions of the form 

_2 

2d 



exp(ic*r)^^ = lira exp ( - ^ c^MT^^z) , (17) 



where c is again some arbitrary vector which does not explicitly depend on R. We can 
always go back to the original monomer coordinates since f = P^^R. The key problem here 
is to find the inverse of M.{^z) in the limit z = {e/a)"^ which depends on all the crosslink 
coordinates {ie,je}eLi- This is accomplished by invoking an identity due to Frobenius, Schur 
and Woodbury [0,113 • An alternative but shorter proof of the theorem is given in appendix 
B. The important finding is that 

M-\z) = -(x-ViV'vy'V') . (18) 

Making use of the fact that = (V^'P)^^V^ is a pseudoinverse of V, i.e. V^V = I, it 
is trivial to show that VV^ is a projector for V. It is also possible to prove the following 
alternative representations \W 

M 
e=l 

Here X = (xi, x^y^) is any orthonormal vector basis Xe, e = 1, M for the M dimensional 
vector space spanned by the pe's in (15). Combining equations (17-19) our central result of 
this section can now be summarized by the following formulas 

exp(zc*r))^ = exp ( - |i c\l - XX') c) (20) 



2 M 

,2 \2 



2d 



c x„ 



e=l 
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It is remarkable that equation (20) is of very simple and special form although the crosslink 
coordinates C = {ie,je}eLi are completely random. X— AfA"* is orthogonal to V and idem- 
potent which can be seen by inspection. As a consequence the only eigenvalues of X — XX'^ 
are and 1 with degeneracies M and — M, and the quadratic form in (20) is semipos- 
itive definite. This assures that the exponent in (20) is never positive. The non-triviality 
of the network problem enters the calculation in form of X which can be found either by 
orthonormalizing the Pe's in (15) or directly from (18). The former is usually accomplished 
by Gram-Schmidt orthonormalization ||13| or numerically by singular value decomposition 
|T6| . Note that equation (18) only requires inversion of an M x M matrix, whereas A4~^ is 
dimensional with N ^ M for a real network. 

It is worthwhile to mention that equation (20) is also the partition function of a Gaussian 
structure (random or not depending on the choice of crosslinks {ig, Jejeli) the presence of 
random external fields bi, ...,hN. This is easily seen by making the transformation in (20) 
c = I?*b back to the original monomer coordinates Rj. In the following section we consider 
more applications of equation (20) with emphasis on RCGNs. 



IV. NON-INTERACTING RANDOM NETWORKS 



The reference quantity in the theory of linear polymers and polymer melts [|11||12| is the 
static structure function of a single non-interacting polymer chain. For a RCGN it is given 
by 

N 

So{k, C) = ( ^ exp (^k(R, - R,)))^ . (21) 

i,j=0 

The structure function can be measured directly in polymer solutions under 0-conditions 
via neutron scattering. Physically it is the Fourier transform of the pair correlation function. 
From equation (20) it is easy to derive an exact expression for 5*0 (k, C) by setting 

c = k(0,. ..,0,1,1,..., 1,1,0,. ..,0)*. (22) 
« ' 

i + 1 to j 
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From equations (20-22) we find 

N 1,2 2 

5o(k,C)= Eexp(-^(|z-j|-(y.-y,)^)). (23) 

i,j=0 

The M dimensional vectors {i = 0, N) are given in terms of tlie ortlionormal basis X 
in (19) 

iy,,...,y^) = iVXy , (24) 

wliere T> is tlie lower triangular matrix defined in (12), and yo = 0. Each y^ is a vector whose 
components depend on the whole set of crosslink coordinates C via X. Equation (23) is of 



similar structure as the result in reference 17 for a RCGN under external stress obtained 



by completely different means. Derivations of structure factors that are based on the affine 



deformation hypothesis can be found in references |T7| , p!8[| but will not be dealt with in 
this work. 

Although our primary interest here are RCGNs, we can apply equation (23) to prob- 
lems with non-random connectivities as well. As a simple example which can be solved 
analytically and to illustrate our formalism, we consider a flexible ring polymer as a triv- 
ial example of a non-random network with only one crosslink. For a closed loop the 
crosslink connects monomer Rq with Rjy, and the crosslink matrix V is of the simple form 
p = pi = (1, 1, . . . , 1)*. Thus X = ^1 = 1/y/N (1, 1, ... , 1)* and with the definition in (24), 
Yi = i/ VN ioT i = 0, N. From equation (23) we get the exact result 

5o(k, Ring) = E exp f - ^-4^(1 - K^)) ■ (25) 



There are more cases in which X can be obtained analytically. Among these are star or 
branched polymers and networks or manifolds with regular, non-random connectivity. More 
details of these applications will be given elsewhere. Since equation (23) is exact, it also 
reduces for M = to the well known Debye function for linear polymers W^ . 



We now specify our discussion to networks with quenched random connectivity. For 
an arbitrary set of crosslinks we are not able to derive analytic expressions for yj in (23). 
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Moreover, within the framework of rephca theory macroscopic physical observables hke 5*0 
are to be averaged over the "frozen" variables, here all possible sets of crosslink coordinates 
C. To make further progress it is therefore absolutely crucial to understand that in the 
thermodynamic limit when and M are sufficiently large any specific but random crosslink 
configuration C would produce the same result for 5*0 if, for example, equation (23) could be 
evaluated analytically. In the literature this is well known as self-averaging. In fact for any 
self- averaging quantity like the free energy, structure factor or radius of gyration (these are 
the quantities we are mainly interested in) performing the quenched average at the end of 
the calculation is not an absolute necessity. To obtain Sq for RCGNs it is therefore sufficient 
to pick one specific set of random numbers < ie.je N, e = 1,...,M, from a suitable 
distribution function and evaluate (23) for this specific but random connectivity. With the 
results from section 3 it is indeed very easy to obtain and carry out the summation in 
(23) numerically (figure 2). To create a specific set of crosslink coordinates C = {ie, je}fli 
we have chosen an uniform distribution P(C) = (1/A^)^*^ for simplicity [Q. Physically 
this means that the frozen crosslinks can in principle be anywhere along the chain without 
restrictions. 

Self- averaging is demonstrated in figure 2 for two random networks with different 
crosslink configurations C and also different network sizes N and M. We find that self- 
averaging is almost perfectly fulfilled even for relatively small networks with ~ 5000 and 
crosslink densities of a few percent. Small deviations between the two networks in figure 
2 are due to finite size effects. The structure function depicted in figure 2 is a universal 
function in a sense that any other RCGN with ~ 5000 led to identical curves. As one 
would have expected Sq does only depend on the number of crosslinks M and monomers 
A^ and not on microscopic details of crosslink positions. In fact our exact calculation shows 
that for RCGNs, S'o(k, C = {ie,je}eLi) = 'S'o(k, c) and that it depends only on crosslink 
concentration c = M/N. The important consequence is that for RCGNs the structure func- 
tion 5*0 (k, c) can be viewed as the equivalent to the Debye function in the theory of linear 
polymers. By construction it is this complete analogy between the theory of linear polymers 
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and our model that will allow us to approach the excluded volume problem in the next 
section very similar as for linear polymers. 

As a by-product of equation (23) we obtain an exact expression for the radius of gyration 
i?g of a RCGN. It is easy to show |1T1 that in the limit k ^ 0, 5'o(k, c)/N'^ = 1 - k^Rl/d, 
and from our numerical calculation we find that = 0.26a^/c, where c = M/N. The 
latter result is also universal in a sense that the prefactor 0.26 ±0.01 ~ 1// was found for all 
networks with functionality / = 4 which were modelled from a uniform crosslink distribution. 
Small fluctuations were mainly due to the finite size of the networks under investigation. 
Similar results have been suggested throughout the literature. However, these calculations 
had to rely upon various approximation schemes. For very large k values of the order 1/a, 
5*0 (k, c) decays as 1/k^ as expected for a non-interacting system. In the intermediate k range 
no simple power-law decay could be found. From the semipositive definiteness of (20) and 
(23) it follows that the scattering function of any simply-connected Gaussian structure is a 
monotonously decreasing function in k^ with a maximum at k = 0. So far further analytic 
progress in (23) depends on simplifying assumptions or approximations about the form of 
Yi's in (24). 

As another illustration of (20) we consider a generalized Gaussian structure in a ho- 
mogeneous electric field E when each monomer is randomly charged with qi = ±q. In 
this case the electric energy is given by — X^i^jERj, and we can directly get the partition 
function from (20). Performing the quenched average over which is analytically possible 
the free energy becomes F(E) = F{0) — {a? q^Y? / 2dkjiT) Tr3^, where Tr means trace, and 
y = V{1-XX^)V\ Applying the definition in (24) we find that Tr 3^ = N{N-l)/2-Y.f=i y- 
and from a numerical calculation TrJ^ ^ 0.5N/c. 

V. RANDOM NETWORKS WITH EXCLUDED VOLUME 
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A. Free energy 



We want to outline how the free energy of a RCGN with excluded volume interaction 
can be obtained without resorting to replica methods. For the above considered network the 
partition function (8) is now expanded in terms of density fluctuations up to second 
order. This approximation is known to be valid in dense systems when fluctuations about 



the mean density are small [rT|. In the Gaussian approximation |n] it is straightforward to 
obtain the free energy density F{v,c) of the interacting system in terms of 5'o(k, c) of the 
non-interacting system 

F{v,c)-F{0,c) = h,T(vp' + J (0?log(l + f^o(k,c))) , (26) 

where p = N/V is the mean physical density. Note that no further averaging over the frozen 
crosslink positions is needed, since 5*0 (k,c) in (23) was already shown to be self- averaging. 
For further reference we quote also the result of the mean field free energy in terms of mean 



density fluctuations ($k)- Up to an irrelevant constant it has been shown [|TT| that 

F{v,c) = \kT J ^5-i(k,c)($k)($-0 , (27) 

where S'~^(k, c) = (v/V) + S'5'^(k, c) is the inverse structure function of the interacting 
system in Gaussian approximation. Higher order terms in the above Landau type expansion 
are readily calculated from equation (20) by modifying c in (22) accordingly, but will be 
neglected in the following discussion. 



B. Crosslinked polymer blends 

Interpenetrating polymer networks (IPNs), semi-IPNs and crosslinked polymer blends 
have broad range of applications in polymer research and material sciences ||T9|j20[| . There 
has been some effort to model these systems using a propagator originally proposed by de 
Gennes which led to physically reasonable results and agreed widely with experiment 
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0]. However, these were semi-phenomenological models and a microscopic theory for multi- 
component networks is to the best of our knowledge still lacking. 

We consider the case of a highly crosslinked two-component polymer blend. For suffi- 
ciently strong incompatibility between the network components we expect phase separation 
on a mesoscopic length scale in which A rich and B rich microdomains are formed. To model 
such a situation we start from a long A-B diblock copolymer with the following structure 

{A monomer if < z < SN , 
- (28) 
B monomer if (j)N < i < N . 

Here and 1 — are the volume fractions of the A and B components. Crosslinks between 
monomers are introduced in the same fashion as in equation (4). The additional crosslink 
between the A and B chain has of course no severe consequences. The interaction is described 
by a proper generalization of Hi in (2) 

I (f>N-l N 

i,j=0 i,j=(j>N 
4>N-l N 

+ 2yAB E E 5(R.-R.i)) • (29) 

We follow here closely Leibler's derivation |]I2[ for diblock copolymer melts. Using the more 
general expression (29) in the Hamiltonian (5), the free energy is calculated in terms of 
one component density fluctuations 5^ = Ht=o~^ exp(ikRj) up to second order. From this 
standard calculation a free energy functional of the same form as (27) is obtained with a 
modifled expression for the structure function S [|l2l . For an incompressible system of sym- 
metric copolymers (f aa = "Web = v and f ab = v + Av) S is given in terms of crosscorrelation 
functions of single component density fluctuations 5'Q^(k, c) = ($k'^^k)o {x,y = A,B) and 
the Flory parameter Xf = pA^; [|12| 

CAA I CBB I O OAB 

S^^S^"" - (50^^)2 N ' ^ ^ 

With the method described in section 4 it is straightforward to calculate S^^ik, c) by a proper 
generalization of (23) for arbitrary volume fraction < < 1. In the disordered phase it 
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is sufficient to consider quadratic terms of the order parameter ($^) in the Landau free 
energy, equation (27). A complete understanding of microphase separation would require 
the study of the homogeneous phase as well as the "ordered" mesophase. As precursor for 
phase instability we analyze the behaviour of the structure function S, i.e., the correlation 
function of the local order parameter, at its maximum. From (30) it is clear that the position 
of the maximum can only depend on the magnitude of k and is independent of Xf- From 
our calculation we find a maximum of S at finite wavevector ko (figure 3). Divergency of 
the structure function in (30) (i.e., 0) at its peak value ko indicates instability of the 

disordered phase and serves as definition for critical Xc (figure 3). More calculational details 



will be presented elsewhere. For a lucid discussion of the above method, see reference [12 



VI. OUTLOOK AND CONCLUSION 

There are many directions to extend this work and open problems in connection with 
macromolecular networks. Here we mention some of them mainly for completeness. 

In this paper we have considered networks that were modelled from a single chain of 
macroscopic size. While this simplification is physically reasonable for highly crosslinked 
polymers ll|-^, it completely fails to describe the correct physics of dilute systems. The 
problem of weakly crosslinked polymers is of special relevance in connection with the vul- 
canization transition, i.e., the transition from a liquid to an amorphous solid state. Upon 
increasing crosslink concentration the gel-fraction of the network grows until at the percola- 
tion threshold an infinitely large network cluster is formed. In this case it was shown that a 
finite fraction of polymers spontaneously localizes with respect to some reference frame 
(e.g., the center of mass of the gel) and therefore can support applied stress. Although we 
are not dealing here with the vulcanization problem (a// monomers are a-priori localized via 
the harmonic potentials, equation (4)), it is instructive to see how the same issue arises in 
our formulation. 

In a dilute network consisting of many unconnected or partially connected polymer 
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chains, each network cluster has its own center of mass mode. In this case it is no longer 
possible to integrate over each center of mass coordinate separately as was done in the be- 
ginning of section 3. The simplifying feature there was of course simply-connectedness of 
the object. Since any further progress required the matrix in (14) to be positive definite, a 
randomly crosslinked dilute network consisting of many chains poses additional mathemati- 
cal problems. It is easy to see that for multiple chains, equation (Al) has to be replaced by 



the more general expression for a polymer melt 



T 



■■■ 
1 ■■■ 



(31) 



^00... 1^ 



Here X„ is the n dimensional unit matrix, n the number of polymer chains, and ® denotes 
a direct product. Equation (31) has a simple physical interpretation. The I's represent the 
backbone of the chain and by removing one element from the diagonal, the chain is cut into 
two pieces, and so forth. It is interesting to note that even for this more general situation, 
it is still possible to "invert" M. analytically, although it is no longer positive definite. This 
was proven in reference by making use of the concept of pseudoinverse matrices which 
leads to a similar result as equation (18). However, the mathematics is more involved and 
will not be presented in this paper. In the framework of replica field theory the vulcanization 
problem was treated in references |^ by invoking a variational ansatz for the localization 
length of monomers. 

Finally we want to make a short comment about the role of entanglements in poly- 
meric networks. Up to now the systems considered were phantom in a sense that the only 
topological restrictions on chain motion were permanent crosslinking junctions modelled 
by (4). In any realistic vulcanization process upon network formation a certain number 
of entanglements is permanently trapped which can be viewed as another form of frozen 
constraint leading to reduction of entropy. However, up to now there are no topological 
invariants known to mathematicians that describe this phenomenon in an unique and rig- 



orous fashion. A simplified version of the entanglement problem was proposed by Edwards 
and coworkers who have modelled entanglements in form of sliplinks P,p3|, i.e., crosslinks 
with the additional freedom to slide along the chain. Formally sliplinks are introduced by 
treating the crosslink coordinates as "hot" variables. To see how equation (3) gets mod- 
ified in the presence of sliplinks, we consider the simplest possible scenario in which all 
crosslinks in (3) are assumed to be sliplinks. To model a more realistic situation we had to 
consider both, crosslinks and sliplinks. Summing (3) over all ie,je independently modifies 
the former crosslink term to {J2i^j=o^ (^i ~ R-j) )*^- By invoking the well known identity 
x^^ = (M\/27ri) § dfj.e'^^ it becomes clear that the sliplink contribution can be 

treated on a similar footing as the excluded volume interaction in (2). A mathematically 
similar but more involved problem arises in the replica formalism, and it had been shown 
lU that the /i-integration is dominated by the steepest descent. Thus in the simple example 
considered the effective excluded volume is reduced by the presence of sliplinks to v — fi. The 
above discussion is also valid for independently distributed crosslinks and sliplinks. Whereas 
the crosslinks modify the Wiener measure to (4), sliplinks renormalize the excluded volume 
parameter. More sophisticated models can be formulated if the degrees of freedom of the 
sliplinks are in addition restricted by the permanent crosslink positions PPB| . 



In summary, we have proposed a microscopic model for RCGNs and its generalization 
to interacting networks and crosslinked polymer blends. As an extension of the Edwards 
Hamiltonian for linear polymers we introduced a similar Hamiltonian for RCGNs taking 
excluded volume interaction into account. We developed a mathematical formalism which 
can solve the non-interacting RCGN exactly. By employing the idea of self-averaging we 
showed how to approach the excluded volume problem in a systematic, perturbative man- 
ner. Our method avoids the well known technical difficulties of replica theory like replica 
symmetry breaking or the n limit. The presented theory provides a new perspective 
on various aspects and open questions of polymer networks with random and non-random 
connectivities that are of general interest. 
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APPENDIX A: 



To prove equation (14) we have to evaluate T>*{zyV+J2^i ^(^e, je))T^ in (9) with matrices 
defined via (10-12). First note that 



^ ••• 0^ 



1 ••• 







(Al) 



Let Bi = (0, 0, 1, 0, 0)* be the ith unitvector of the canonical basis in + 1 dimensional 
space. Then 

M M 

e=l e=l 
M 

= E (e. 

e=l 



e,e)K-Se) 

= {K - e^J, (e^^ - e^.^)) ((e,^ - e^J, (e,^ - e^.^ 
From (A2) we obtain equations (14) and (15) by setting Pe = — 'D*(ej^ — e^-^). 



(A2) 



APPENDIX B: 



We want to find the inverse of the matrix defined in equation (14), 

M{z) =z{x+-VV') , 



(Bl) 



in the limit z — > 0. The difficulty with (B3) is that X cannot be neglected against the 1/z 
term, since the crosslink contribution alone is highly singular with degeneracy N — M. We 
proceed by writing the inverse of A4. quite formally in terms of its Neumann series 
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1 1 — 1 1 °° 1 

M-\z) = -(I+-VV') = - E ( - -'P'pT 

z > z^^ z 

1 / 1 OO 1 



VM . (B2) 



2; \ z ^ z 

With the definition in (15) it is easy to show that for tetrafunctional crosshnking junctions, 
the crosshnk matrix V is of full rank M. It is well known that in this case the normal form 
V^V is positive definite [|13| and that therefore its inverse exists. Using this information we 
obtain equation (18) from (B4) by letting 2; 0. Convergency of the above manipulations 
is proved by multiplying the final result for TW^^ with hA in (B3) which gives identity. This 
also verifies the correctness of the result in (B4). 
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FIGURES 

FIG. 1. A single-chain network. Only the beads at the crosslink sites are shown. For high 
crosslink concentrations far above the percolation threshold a single-chain network is practically 
indistinguishable from a network that is modelled from many polymer chains. 

FIG. 2. Scattering function So{q, C) for two different RCGNs derived from (23). = k'^a'^/2d; 
c = M/N; solid lines N = 10000, c = 0, 1, 2, 3, 4, 5% (left to right); open circles N = 5000 ( 
only shown for c = 2, 4%). C = {ie, 3e}e=i were picked from the interval [0, A^] at random. The 
orthonormalization of V was performed with a standard singular value decomposition algorithm 
[16]. Due to self-averaging curves with same crosslink concentration (here c = 2, 4%) are identical. 
For c — ^ we recover the Debye function for linear polymers. 

FIG. 3. Structure function (30) of a symmetric crosslinked polymer blend for three different 
values of the Flory parameter xf = 0, 2c, 4c (upwards) with c = M/N = 0.02 (M = 200, 
N = 10000) and <p = 0.5. We find that a maximum of S occurs at finite wavevector qo = 0.23 
(g^ = k^o^/2(i). For the above values of and c the disordered phase becomes unstable for 
Xf > Xc ~ 4.3c. 
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